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ANALYTICAL  AND  EXPERIMENTAL  INVESTIGATION  OP  PERFORMANCE 
OP  SUPERSONIC  EJECTOR  NOZZLES 

Shen  Hull!  and  Chen  Zbongqing 

This  paper  presents  an  analytical  method  for  calculating  the 
flow  field  and  performance  of  supersonic  ejector  nozzles.  The  cal¬ 
culations  involve  the  real  sonic  line  at  the  exit  of  the  primary 
nozzle,  the  inviscid  primary  flow  field,  the  correction  for  the  vis¬ 
cosity  effect,  the  pumping  characteristic  and  the  thrust  character¬ 
istic. 

In  order  to  have  the  calculated  results  agree  with  the  experi¬ 
mental  data,  the  real  sonic  line,  instead  of  the  plane  sonic  line, 
is  taken  as  the  initial  base  line  of  calculations.  The  real  sonic 
line  is  obtained  by  joining  the  points  of  intersection  of  constant 
flow  angle  lines  (isoclines)  and  the  throat  region  with  Mach  lines 
at  the  lip  of  the  primary  nozzle. 

The  inviscid  primary  flow  field  of  the  nozzle  is  calculated  at 
first  and  then  a  correction  is  made  to  account  for  the  viscous  effect. 
The  method  of  correction  for  the  viscous  effect  proposed  in  this 
paper  replaces  the  original  geometric  coordinates  of  the  ejector 
shroud  with  the  corrected  geometric  coordinates  which  are  obtained 
by  superimposing  on  the  original  geometric  coordinates  the  displace¬ 
ments  of  the  mixed  region  and  the  boundary  layer.  On  the  basis  of 
the  corrected  coordinates,  the  actual  primary  flow  field  and  pumping 
performance  of  the  nozzle  are  then  calculated.  The  proposed  method 
proves  to  be  quite  simple  and  accurate  . 

Calculations  were  performed  on  a  320  digital  computer  and  model 
tests  were  conducted  on  a  ground  test  facility.  The  analytical  and 
experimental  results  are  found  to  be  in  fairly  satisfactory  agreement. 

For  a  relatively  high  secondary  flow  rate,  the  primary  flow  con¬ 
tinues  to  accelerate  and  expand  in  the  divergent  section  of  the  "fluid 
wall  boundary"  after  leaving  the  primary  nozzle,  while  the  secondary 
flow  accelerates  in  the  passage  formed  between  the  shroud  and  the 
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FLOW  CONFIGURATION  [1]  [2] 


Ejector  nozzle  with  high  secondary  flow  rate  (Figure  1) 


Figure  1.  Flow  configuration  in  ejector  nozzle  with  high 

secondary  flow  rate 

Key:  1 — secondary;  2--primary 

"fluid  boundary".  Depending  on  the  dimension  of  the  nozzle  and  the 
upstream  flow  conditions,  there  are  three  possibilities: 
the  secondary  flow  remains  subsonic; 

the  secondary  flow  becomes  sonic  at  the  exit  plane  of  the  nozzle; 
the  secondary  flow  becomes  supersonic  at  the  exit  plane  of  the  nozzle. 

Due  to  the  viscosity  of  the  gas,  a  mixing  layer  of  variable 
pressure  is  formed  between  the  primary  flow  and  the  secondary  flow 
and  a  boundary  layer  is  formed  along  the  shroud.  However,  the  vis¬ 
cous  effect  is  relatively  small  for  a  high  secondary  flow  rate. 

For  a  very  low  secondary  flow  rate  (approaching  zero),  the  pri¬ 
mary  flow  practically  expands  as  a  free  jet  into  a  region  of  con¬ 
stant  pressure.  The  free  Jet  boundary  almost  sticks  to  the  divergent 
section  of  the  shroud  and  is  recompressed  to  produce  a  shock.  Up¬ 
stream  of  the  compression  region,  the  low  speed,  constant  pressure 
flow  restricted  by  the  constant  pressure  Jet  boundary  forms  a  dead 
region.  Due  to  the  viscous  effect,  a  mixing  region  forms  along  the 
boundary  of  the  primary  flow.  Mass  and  energy  transfers  between  the 


Flow  configuration  with  secondary  flow  rate  (Figure  2) 


Figure  2.  Flow  configuration  of  the  ejector  nozzle 
with  low  secondary  flow  rate 
Key:  1 — secondary;  2 — primary 

primary  and  the  secondary  flows  occur  in  this  mixing  region.  The 
primary  flow  continuously  drags  out  a  certain  flow  rate  from  the 
dead  region.  This  flow  rate  is  provided  by  the  supply  system  of 
the  secondary  flow. 


METHOD  OF  CALCULATION 


1.  Calculation  of  the  real  sonic  line  C3] 


The  real  sonic  line  must  be  determined  after  the  flow  angles 
at  the  transonic  throat  region  of  the  primary  nozzle  and  the  Mach 
lines  generating  from  the  lip  of  the  primary  nozzle  are  calculated 
(Figure  3).  Assuming  that  the  flow  is  steady,  irrotational,  isen- 
tropic  and  two-dimensional,  it  can  be  described  in  terms  of  the 
stream  function  by: 


(  »’ }  *  •*’  ]■•[*(*>')  ] 9n 


<  i  > 


where  the  subscripts  denote  the  derivatives  in  the  x  and  y  direct¬ 
ions.  The  following  velocity  variables  are  introduced: 

*-<•»♦•*)*  <»> 

<*> 

Then  equation  (1)  becomes: 


3 


U’j y.' -$)&••■ 

Again  introducing  the  transformed  velocity  defined  by: 

if  «»' 

the  above  equation  is  simplified  into 


where 


KC-j-O/n  -w> 


Upon  introducing  the  "tangent  gas  approximation",  i.e., 

K(M)  *  1,  the  hodograph  relation  (6)  is  simplified  into  the  Cauchy- 
Riemann  equation.  Hence,  the  flow  field  can  be  solved  by  complex 
variable  techniques.  The  incompressible  and  compressible  flows  are 
related  by  the  transformation 


upon  direct  integration 


1  (  )*' 

.•Lizr* 

"  1  *  l  -  Ml 


where  q  and  M  represent  the  velocity  and  the  Mach  number  at  the 
m  m 

matching  state.  Since  the  Cauchy-Riemann  condition  is  satified,  a 
set  of  sources  and  sinks  satisfying  the  required  boundary  conditions 
is  employed  to  determine  the  complex  potential  function: 

•/*) 

Ftm- .#)-(.  -  *  ■  <»> 

where  .  Upon  differentiating  the  above  equation,  the  com¬ 

plex  velocity  is  obtained: 


where  ./«*•*.  • 
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Upon  applying  this  solution  to  the  conical  nozzle  and  intro¬ 
ducing  the  complex  variable  z«x+iy,  the  hodograph  plane  can  be 
transformed  into  the  physical  plane.  Hence 


4* ' 


4F 


where  the  bar  over  the  last  term  denotes  that  the  complex  conjugate 
is  to  be  taken.  Integration  of  equation  (11)  gives  the  physical 
location  (x,y)  of  any  point  In  the  hodograph  plane.  The  constant 
of  integration  can  be  determined  by  assigning  the  physical  location 
(0,1)  to  correspond  to  the  nozzle  lip  (w^,  a)  in  the  hodograph  plane. 
The  distribution  of  the  flow  angle  in  the  throat  region  is  particular¬ 
ly  useful.  The  lines  of  constant  flow  angle  are  also  called  "iso¬ 
clines".  In  order  to  obtain  the  equation  for  the  isoclines,  equa¬ 
tion  (11)  is  Integrated  numerically  from  w_  to  u),  for  selected  values 

a  j 

of  the  flow  angle  between  0  and  a. 


4P 


(ID 


The  method  of  characteristics  is  employed  to  calculate  the  Mach 
lines  generated  from  the  nozzle  lip.  The  characteristic  equations 
are: 


J  p  «»n#  •  iv 

^  y 

•It  '»  >.  J  *  -  \r  i>  U 


( i.J) 


where  y  is  the  local  Mach  angle.  The  static  pressure  p  and  the  flow 
angle  9  are  chosen  as  the  independent  variables.  The  expansion  at 
the  nozzle  lip  is  of  the  centered  wave  type  and  satisfies  the 
Prandtl-Meyer  function.  The  flow  at  the  nozzle  lip  turns  through 
the  Prandtl-Meyer  angle,  and  this  expansion  process  is  divided  into 
a  finite  number  of  discrete  steps.  The  first  point  on  the  sonic 
line  is  the  point  of  intersection  of  the  first  isocline  and  the  first 
Mach  line  originating  from  the  nozzle  lip  (Figure  3).  The  locations 
of  the  other  points  on  the  sonic  line  are  determined  In  the  same 
fashion. 

2 .  Calculation  of  the  lnylscld  flow  field  and  the 
characteristics  for  high  secondary  flow  rate 
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Figure  3.  Real  sonic  line 

Key:  1 — jet  boundary;  2 — isoclines;  3 — sonic  line 

(1)  Assumptions : 

The  secondary  flow  is  a  one-dimensional  isentropic  flow;  the 
static  pressure  along  the  "fluid  boundary"  between  the  primary  flow 
and  the  secondary  flow  is  common  on  both  sides;  the  primary  flow  is 
inviscid  and  is  subjected  to  the  influence  of  the  characteristics. 

(2)  Given  conditions: 

stagnation  pressure,  stagnation  temperature  and  mass  flow  rate  of  the 
primary  flow; 

stagnation  pressure  and  stagnation  temperature  of  the  secondary  flow; 
ambient  pressure; 

geometric  dimensions  of  the' ejector  nozzle. 


(3)  Principal  equations  applied  [5]: 


Figure  4.  Calculation  scheme  for  the  ejector  nozzle 
Key:  1 — characteristics;  2 — determined  points;  3 — streamline; 
4 — calculation  point;  5— axis  of  revolution 
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(4)  Calculation  procedures:  An  arbitrary  secondary  mass  flow 
rate  Is  taken.  Based  on  the  condition  of  equal  static  pressure  for 
the  primary  flow  and  the  secondary  flow  at  the  exit  of  the  primary 
nozzle,  the  expansion  fan  is  determined.  The  location  of  the  "fluid 
boundary"  is  determined  by  iteration.  The  continuity  equation  is 
written  for  each  section  along  the  secondary  flow,  while  the  equa¬ 
tion  for  the  characteristic  originating  from  this  point  is  written 
at  the  same  time.  The  calculation  is  conducted  until  the  nozzle 
exit  is  reached.  If  the  static  pressure  of  the  secondary  flow  at 
the  nozzle  exit  is  found  to  be  the  same  as  the  ambient  pressure,  the 
secondary  flow  is  subsonic.  If  In  case  the  secondary  flow  is  found 
to  be  supersonic  at  the  exit,  then  there  exists  a  throat  section  In 
the  secondary  flow  where  the  Mach  number  is  1.  The  mass  flow  rate 
selected  must  then  be  modified  such  that  it  is  equal  to  the  choked 
flow  rate . 

3.  Correction  for  viscous  effect  [4],  [3],  [7] 

Reference  [4]  included  the  viscous  effect  by  correcting  the 
cross  sectional  area  of  the  throat  of  the  secondary  flow.  [3],  on 
the  other  hand,  corrected  the  boundary  between  the  two  flows  and 
included  the  boundary  layer  along  the  shroud.  The  former  is  simpler 
but  less  accurate,  and  the  latter  method  is  more  complicated.  The 
present  paper  proposes  the  shroud  coordinates  correction  method, 
which  can  be  described  as: 

(1)  calculation  of  the  inviscid  flow  field  with  the  initial  dimen¬ 
sion  of  the  shroud  of  the  ejector; 

(2)  calculation  of  the  displacement  thickness  of  the  mixing  region 
between  the  primary  flow  and  the  secondary  flow; 
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(3)  calculation  of  the  boundary  layer  displacement  thickness  along 
the  inner  wall  of  the  shroud; 

(4)  correction  of  the  throat  section  of  the  secondary  flow  by  sub¬ 
tracting  the  boundary  layer  displacement  thickness  from  the  orig¬ 
inal  throat  section  dimension  and  adding  to  it  the  mixing  layer  dis¬ 
placement  thickness; 

(5)  based  on  the  corrected  shroud  coordinates,  the  entire  flow  field 
is  calculated  again  by  the  inviscid  method.  The  analysis  is  des¬ 
cribed  in  detail  below. 


Figure  5-  Velocity 
profile  in  the  mix¬ 
ing  region 


.  .  # 

(1)  Calculation  of  the  displacement  thickness  6-.  of  the 

secondary  flow  in  the  mixing  layer  1 

If  the  mixing  layer  is  assumed  to  be  a  constant-pressure  region, 
the  velocity  profile  in  the  mixing  layer  (Figure  5)  can  be  obtained 
by  the  following  equation  [4]: 

* .  I - '  1  ■/*  ( y  v**  •  w)  - 1 7*  * '  7**rl  •  _ 

where 

u^  -  velocity  of  the  secondary  flow 
u  -  velocity  of  the  primary  flow; 
u  -  velocity  component  in  the 

x  direction  «<’*  *  \  f 

«f  *  *  ^ 

i  * '7 '  a  -  similarity  parameter,  given  by  a  semi-empirical 

formula;  y,x  -  physical  coordinates  in  the  mixing  layer; 
velocity  at  the  point  y  =  0:  u  (that  is  *  *' )' 

n  is  the  similarity  coordinate,  that  is,  for  equal  values  of 
then  the  jj  value  is  also  the  same.  In  other  words,  the  velocity 
profile  ataeach  section  of  the  mixing  layer  is  similar. 
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Upon  taking  the  Prandtl  number  as  1,  neglecting 
assuming  C  to  be  constant,  the  temperature  profile  can  be  obtained 
from  the  velocity  profile  directly,  that  is 


di 


and 


7  • 

'  ;  7: 


w.[(£  '•)•('-  fit* : 


‘il> 


*  * 

where  T  -  stagnation  temperature  of  the  mixing  region;  T  -  stag- 

*  a 

nation  temprature  of  the  primary  flow;  -  stagnation  temperature 
of  the  secondary  flow. 


Given  the  velocity,  stagnation  temperature  and  static  pressure 
(equal  to  the  static  pressure  of  the  primary  and  secondary  flows  at 
the  boundary  of  the  mixing  layer)  at  each  point  of  the  mixing  layer, 
the  displacement  thickness  of  the  mixing  layer  can  be  obtained. 

The  velocity  profile  and  the  temperature  profile  above  are  given 
in  terms  of  the  internal,  coordinate  system  (x,y)  which  is  related  to 
the  reference  coordinate  system  (X,Y)  by 

Y  =  y-y 

m 

The  jet  boundary  in  the  mixing  layer  is  denoted  by  j  (see  Figure 
5).  The  points  for  which  y=y^  separate  the  mixing  layer  into  two 
regions:  the  mass  flow  rate  in  the  range  from  n j  to  nRa  is  exactly 

the  portion  of  the  mass  flow  rate  of  the  primary  flow  entering  the 
mixing  layer;  and  the  mass  flow  rate  between  q,  and  nRb  is  exactly 
the  portion  of  the  mass  flow  rate  of  the  secondary  flow  entering  the 
mixing  layer. 

•>,  is  solved  by  the  principle  of  conservation  of  mass  and  momen¬ 
tum  before  and  after  the  mixing: 

f  ’  f  .<•/,>  ?  t  i  ( 2  'J. ) 


where  c  =  Crocco’s  number. 
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ym  is  also  determined  by  the  principle  of  conservation  of  mass  and 
momentum.  In  terms  of  the  similarity  coordinate: 

-  -j  iU> 


Analogous  to  the  displacement  thickness  concept  of  the  boundary 

« 

layer,  the  secondary  flow  displacement  thickness  6  ^  of  the  mixing 
layer  is  found  to  be 


where 


v  a  1  — 

n 

2.758  -C., 


■  I2> 


v  <  i  -c:,r¥^  ’ 


O 

o7 


jff' 


i 


l-c. 


r  «  CU  1-r!) 


cn> 

<2t) 

<2*> 

(30) 


Key:  1 — mix;  2 — secondary;  3 — primary 


£ 

(2)  Calculation  of  the  displacement  thickness  6.  ,  of  the 
boundary  layer  along  the  inner  wall  of  the  shfroud  D *i* 


£ 

The  effect  of  the  boundary  layer  displacment  thickness  6.  n 
along  the  inner  wall  of  the  shroud  on  the  secondary  flow  rate  is 
estimated.  The  properties  of  the  boundary  layer  are  calculated 
based  on  the  equivalent  length  which  is  defined  as: 


whei»e 


}  • 


u  y 
l  *0.2 VV 


(311 


$  is  taken  as  1.20  or  1.25,  depending  on  the  Reynold's  number. 
Given  the  stagnation  pressure  and  the  stagnatiOi  temperature,  the 
local  Reynold's  number  can  be  calculated: 


where 


•  —  ^f7-  MU  1 

C«> 

. .  ....  -<r*> 

•9.5./  R  — 

(39> 

» 

The  displacement  thickness  5^  ^  and  the  momentum  thickness  0  can 

then  be  calculated  in  terms  of  the  local  Reynold's  number  R  .  If 

6  x 
Rx  is  of  the  order  10°,  they  are  found  to  be 

*•<  i  *0.8 Ml >*"**•  w*’**  <*** 

»  .0.034-- •<  I  ♦O.l  <*> 

1— b.l. 
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If  Rx  is  of  the  order  10  ,  they  are. expressed  by: 

O’^O.OOt-  *-<10  +  0. *.!/!)*■•  •  <.1*> 

«  ^0.022- *•<!. 0  +  0. U/J> — '  UD 


1— b .  1 . 


The  above  equations  are  valid  for  K  =  1.4. 


Since  the  static  pressure  in  the  boundary  layer  is  assumed  to 

* 

be  a  constant,  the  pressure  on  the  surface  6  is  always  taken  as  a 
criterion  for  iteration  during  the  calculation.  Hence,  the  calcula¬ 
tion  for  the  points  inside  the  boundary  layer  is  of  the  same  stand¬ 
ard  form  as  that  for  the  points  on  the  surface. 


4 .  Calculation  of  the  thrust  parameter  [5] 


The  thrust  of  the  engine  is  expressed  by  the  followine  relation 


Gg  f.  .  p,A.  -  p,  *.-y.  •#.-  gV< 


where  R  is  the  total  thrust 
a 

*  f,  ■  p.  •  a. 


•  m 


According  to  the  momentum  theorem,  the  above  equation  can  be 


rewritten  as 


in  which 


-F  •  F*  -  -  Ph  A% 

<  :o> 

/  ,H  1 .0  *•  K .Cm  *  C » •  .V!  M.  •  /•ttlttt'f.  -1) 

<  tl  > 

r. - p. .  .Mi.!!* k.‘  wr  >c» 

<  12; 

•V.«  -|"«  v  ’•>!  •• 

]  it 

}  % 

where  C 
C 
C 

y 


f 

p 

v 


fictitious  coefficient  of  the  shroud 
flow  rate  coefficient  of  the  primary  nozzle 
velocity  coefficient  of  the  primary  nozzle 
radius  of  the  shroud  wall 
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If  the  viscous  effect  is  neglected,  the  flow  rate  coefficient 
and  the  velocity  coefficient  of  the  primary  nozzle  are  obtained  by 
the  following  expressions  respectively: 

f:*X  -mAJy-  tia.Wt) 

'•  ii*yJy 

)  i*y<*m4dy  -  nn4J*) 

in  which  the  integrations  are  taken  along  the  sonic  line. 

If  the  viscous  effect  is  taken  in  account,  the  flow  rate  coeff¬ 
icient  and  the  velocity  coefficient  of  the  primary  nozzle  are  found 
by: 

r.  '  “*> 


r;  i  \  (  *  in /*«/•*)  u  ... 

#*  *  K  .•!'  ’**)  j  >  -  yf  «>*  ii.si-iJr)  * 


where 


K«  -o  m 


and  <j>c  -  semi-vertex  angle  of  the  primary  nozzle 

<t>  -  angle  between  the  flow  velocity  on  the  sonic  line 

and  the  x-axis 

y  -  radial  coordinate  of  the  points  on  the  sonic  line 
x  -  axial  coordinate  of  the  points  on  the  sonic  line 
Rc  -  Reynold's  number 

COMPARISON  BETWEEN  THE  NUMERICAL  AND  THE  EXPERIMENTAL  RESULTS 

A  program  for  the  numerical  method  proposed  here  has  been  devel¬ 
oped  and  the  flow  field,  pressure  distribution  along  the  nozzle  wall 
and  the  pumping  characteristic  of  the  ejector  nozzle  are  calculated. 
The  comparison  between  computed  results  and  experimental  data  are 
summarized  below: 

1.  The  flow  field  and  the  pressure  distribution  along  the  wall 
of  the  ejector  nozzle 


The  input  data  are: 

18  total  number  of  node  points  on  the  shroud 

20  number  of  points  taken  along  the  sonic  line 

29  number  of  expansion  waves  taken  as  the  characteristics  net 
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N4  =  100  upper  limit  of  the  iterative  variable  J 
XW  *  0  x  coordinate  of  the  shroud  at  the  exit  plane  of  the 
primary  nozzle 

XC  =  0  x  coordinate  of  free  jet  at  the  exit  plane  of  the 
primary  nozzle 

YC  =  0.045  m  radius  of  the  primary  nozzle  at  the  exit  plane 
V  =  0.88  kg/s  mass  flow  rate  of  the  secondary  flow 
G  =14.4258  kg/s  mass  flow  rate  of  the  primary  flow 

CL 

p£  =  30702  kg/m  stagnation  pressure  of  the  secondary  flow 
T*  =  288  K  stagnation  temperature  of  the  secondary  flow 
=  1.4  specific  heats  ratio  of  the  secondary  flow 
=  29.27  kg-m/kg  K  gas  constant  of  the  secondary  flow 
T*  =  288  K  stagnation  temperature  of  the  primary  flow 

cz 

p*  =  97992  kg/m  stagnation  pressure  of  the  primary  flow 
R  =  29.27  kg-m/kg  K  gas  constant  of  the  primary  flow 

cL 

X^  =  0.090  m  distance  from  the  exit  of  the  primary  nozzle  to  the 
exit  of  the  shroud 

K  =  1.4  specific  heat  ratio  of  the  primary  flow 

cL 

Pg  =  10333  kg/m  ambient  pressure 
Input  data  for  the  shroud  profile  are: 


Y  >  »  X/YC 

&3*Y/YC 

xs*x/rc 

0 

1.2780 

0.6889 

1.12S9 

ft.  111! 

1.2438 

0.^516 

1.1334 

0.2272 

.  2133 

0.8222 

1.1421 

U.3333 

1.1809 

0.8889 

1.1542 

0.3778 

1.1680 

1.1111 

1.1933 

0  4444 

1.1311 

1.2324 

‘i.  m* 

I.I4IV 

1 . 355f 

1.2710 

••  >'56 

1.1341 

i.TW 

1.310* 

..  /  »<»<) 

1.1263 

1.35CU 

The  output  data  are  the  pressure  p,  the  coordinates  x,y,  the 
flow  angle  <J>,  mass  rate  of  flow  of  the  secondary  flow  GS  and  other 
necessary  parameters  at  all  points  of  the  characteristics  net. 

The  characteristics  of  the  ejector  nozzle  obtained  are  shown 
in  Figure  6. 
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Figure  6.  Characteristics  of  the  ejector  nozzle 


/  / 


The  computed  values  and  the  experimental  data  of  the 
pressure  distribution  along  the  shroud  are  shown  in  Figure  7 


Figure  7.  Pressure  distribution  along  the  shroud 
Key:  1 — experimental  data;  2 — calculated  curve 


2 .  Pumping  characteristic  of  the  e.lector  nozzle 

In  Figure  8,  curve  1  is  the  computed  result  for  inviscid  flow, 
and  curve  2  is  the  corrected  result  for  viscous  flow.  The  figure 
indicates  that  the  computed  result  after  correction  for  the  viscous 
effect  is  in  good  agreement  with  the  experimental  data. 
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Figure  8.  Pumping  characteristic  of  the  ejector  nozzle 
Key:  1 — secondary;  2 — primary;  3 — experimental  data;  4 — calculated 
values  of  inviscid  flow;  5 — calculated  values  of  viscous  flow 
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Figure  8  also  indicates  that  for  the  condition  of  high  second¬ 
ary  flow  rate,  the  error  between  the  inviscid  result  and  the  experi¬ 
mental  data  is  no  more  than  10% . 

CONCLUSION 

1.  The  numerical  method  proposed  is  found  to  be  satisfactory. 

2.  The  correction  for  viscous  effect  proposed  is  relatively 
simple  but  still  accurate. 

SYMBOLS 

-  stream  function 
p  -  density 
a  -  speed  of  sound 

u  -  velocity  component  in  the  axial  direction 

v  -  velocity  component  in  the  radial  direction  or  flow  velocity 
at  the  exit  of  the  nozzle 

p  -  pressure 

U  -  Mach  angle 

<|>  -  flow  angle  (angle  between  the  streamline  and  the  horizontal  axis) 

0  -  angle  between  the  characteristic  and  the  radial  axis,  or  the 
momentum  thickness  of  the  boundary  layer 

X  -  equivalent  length 

y  -  radial  distance  normalized  by  the  radius  of  the  primary  ndzzle 
x  -  axial  distance  normalized  by  the  radius  of  the  primary  nozzle 
u)  -  transformed  velocity 

a  -  semi-conical  angle  of  the  primary  nozzle  at  the  exit 
T  -  temperature 
G  -  mass  flow  rate 
M  -  Mach  number 

X  -  velocity  coefficient  of  the  flow 
K  -  specific  heat  ratio 
R  -  thrust  or  gas  constant  or  radius 
A  -  area 

F  -  impuse  or  functional  relation 
n  -  coordinate  along  the  Mach  line 


5  -  coordinate  along  the  Mach  line 

B  -  constant  determined  by  the  Reynold's  number 
Y*  -  coefficient  of  viscosity  expressed  in  terms  of  stagnation 
parameters 

d*  -displacement  thickness  of  the  boundary 

D  •  X  • 

Subscripts : 

primary,  a  -  primary  flaw 
secondary,  b  -  secondary  flow 
C  -  exit  of  the  primary  nozzle 
e  -  exit  plane  of  the  ejector  nozzle 

H  -  ambient  condition 

G  -  total 

C  -  total 

0  -  inlet  of  the  ejector  nozzle 

SH  -  shroud 

Superscript:  *  -  stagnation 
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